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Two atoms in an optical lattice may be made to interact strongly at higher partial waves near 
a Feshbach resonance. These atoms, under appropriate constraints, could be bosonic or fermionic. 
The universal I = 2 energy spectrum for such a system, with a caveat, is presented in this paper, and 
checked with the spectrum obtained by direct numerical integration of the Schrodinger equation. 
The results reported here extend those of Yip for p-wave resonance (Phys. Rev. A 78, 013612 
(2008)), while exploring the limitations of a universal expression for the spectrum for the higher 
partial waves. 
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I. INTRODUCTION 



Normally, the scattering between two neutral atoms in 
higher partial waves is suppressed due to the centrifugal 
barrier. The situation may change, however, by sweeping 
a static magnetic field near a Feshbach resonance (FR) at 
a higher partial wave [H, Ql- In such a situation, the rela- 
tive energy of the scattering atoms is at near-coincidence 
with the energy of a quasi-bound molecular state with 
a nonzero angular momentum. FR was first observed 
with bosons interacting in the s-state 0], and later in 
higher partial waves Two spin polarized fermions 
may interact in a relative odd-^ state in a single channel 
ultra-cold gas, or in an eyen-l state when the two atoms 
arc in distinct spin states. The single-channel p-wavc FR 
was first observed by Regal et al. [l| between "^^K atoms, 
while the fermionic s-wave between these atoms in dis- 
tinct spin states was observed by the same group a little 
earlier [J. 

The two-particle energy spectrum in a spherical har- 
monic oscillator (HO) near a s-wave FR is found to be 
universal when the range of the interaction between the 
two atoms is much smaller than the average interparticle 
distance in the HO, i.e. the oscillator length [a, @, 01- 
Near the Feshbach resonance, this spectrum has been 
checked experimentally Q in the limit of low tunneling 
for fermionic "^^K atoms in an optical trap in two dis- 
tinct spin states. This is remarkable because the theo- 
retical spectrum is obtained in a one-channel approxima- 
tion, and fits the data even with the first term in the 
effective range expansion. More recently, Yip has ob- 
tained the spectrum of two identical fermions in a HO 
near a p-wave FR. It is known that the p-wave FR gets 
split [lOl between \mi\ = 1 and 0. Yip's calculation is 
for mi = 0, but may be generaised for nonzero mi. For 
higher partial waves, however, it was stated that the en- 
ergy spectrum could not be expressed in terms of the 
parameters in the scattering amplitude and the oscilla- 



tor constant 0, [Tl|. In this paper, using a method first 
formulated by Jonsell @ for Z = 0, we express the energy 
spectrum for the higher partial waves in terms of the scat- 
tering parameters. We examine the I = 2 case in detail, 
and compare the spectra obtained from this analytical 
formula with the results obtained by a numerical inte- 
gration of the Schrodinger equation. Yip's assertion is 
justified, but only in a narrow region at resonance where 
the ground state energy goes to zero. The analytical for- 
mula accurately reproduces the excited state spectrum 
in the entire range, but is shown to be not accurate for 
the ground state in a narrow band across the resonance. 
With this exception, even though the energy spectra for 
/ = 1 and I = 2 are extremely sensitive to the choice of 
the effective range parameter, we find that the levels are 
practically unchanged by introducing the next term in 
the expansion, which is shape-dependent. 

In calculating the p wave spectrum, Yip used the ef- 
fective range expansion for higher partial waves [T^ 
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with 1 = 1. Note that the scattering length ai and 
the effective range r/ have the dimensions of (1/)^'+^ 
and (L)~^'+^ respectively. For single channel elastic 
scattering by a power-law potential r~", and without 
any virtual transition to other channels like an excited 
quasi-molecular state, ai in Eq.([T]) is only defined if 
n > {21 + 3) [13. For the effective range r; to exist, 
the restriction is even more severe, n > (2^ + 5). Taking 
n = 6 for the asymptotic behaviour of the interatomic 
potential, we see that for 1 = 1, only ai exists, but not 
ri. For 1 = 2, neither a2, nor r2 is defined [H, 

The situation changes, however, for dressed atoms near 
a Feshbach resonance. Consider scattering in an open 
channel which is coupled to a resonance in the closed 
channel through a spin-dependent two-body interaction 
W that depends on the relative distance r between the 
two atoms [l^. One part is the long range tensor inter- 
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action between the two dipoles that falls off as . The 
other is the spin-exchange interaction between the two 
valence electrons of the alkali atoms that are closer than 
the uncoupling distance r„. For r < r„, the nuclear and 
electronic spins that were otherwise coupled in an iso- 
lated atom get uncoupled due to one atom's proximity 
to the other atom. The resulting spin-exchange interac- 
tion has a range r„ which is smaller than the so-called 
van der Waals distance, a measure of the length scale for 
the potential. The coupling potential W may there- 
fore be regarded as of shorter range than in cases 
where the spin exchange interaction is dominant. 

Elimination of the closed channel results in an addi- 
tional effective potential in the single-channel formalism. 
This effective potential, which dominates the open chan- 
nel scattering near the Feshbach resonance, contains W 
quadratically. The effective range expansion should be 
valid for the higher partial waves for this short range 
potential, as will be discuused in the next section. 



II. THE TWO-BODY PROBLEM 



A. Two-channel scattering 



energy-dependent scattering length that dominates near 
the Feshbach resonance. Since we have in mind partial 
waves / > 2, and Vop{r) = -(fi^ /M)Ce/f , the back- 
ground scattering length does not exist [ij, . Indeed, 
the phase shift caused by this potential is vanishingly 
small for small k. This background phase shift, denoted 



by<5f^(fc), is given by [l5 

tan(5f3(fc) ^ |2-5Ce 



r{5)T{l - 3/2) 
r2(3)r(Z + 7/2) 
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The distortion in the incident wave 0^ for small k may 
therefore be neglected, and the scattering for the higher 
partial waves at low energies is governed by the energy- 
dependent short range potential Veff. As usual, the en- 
ergy denominator in it may be related to the Zeeman 
splitting due to a sweeping magnetic field, giving rise to 
a large variation of the scattering length a; near FR. We 
may then express the flow of the energy levels as a func- 
tion of the scattering length, assuming a fixed effective 
range r; in the expansion ([1]). This we proceed to do, 
generalising a method first proposed by Jonsell Q for 
I = 0. We confine our treatment to mi — 0. 



We first recall the well-known results of two-channel 
scattering, following the treatment and notation of 
Cohen- Tannoudji The Hamiltonian is a 2 x 2 

matrix, with the diagonal elements Hop and Hd, and the 
nondiagonal coupling given by the short range spin ex- 
change potential W. Here Hop = (— -I- Kp), with Vop 
given by the van der Waal power law potential whose 
asymptotic fall-off goes like r~^. The closed channel 
Hamiltonian, in the single-resonance approximation, may 
be expressed as Hd = Eres\4'res >< (f'resl, with the res- 
onance energy at Eres- 

The wave function -0 is a column matrix with two com- 
ponents (f>op and 0c/- It is straightforward to show that 
the two-channel scattering problem may be described en- 
tirely in the open channel by a Lippman-Schwinger equa- 
tion in which the incident wave is distorted by the open 
channel potential Vop, and the scattering is by the effec- 
tive potential Veff. More explicitly, the (relative) out- 
going wave of the dressed atoms {(pop >, in terms of the 
distorted wave|0^ >, is given by 

\(bop >= 10^ > +GtpiE)Vff\<Pt > , (2) 

where 

y ^ ^ \'i'res >< 4>reA ^ 

E - Eres- < (l)res\WG+p{E)W\<j>res > 

(3) 

In the above, G+p{E) = [E - Hop + ie)"^ The distorted 
incident wave > itself obeys a Lippman-Schwinger 
equation in which the incident wave is a plane wave, 
the Green's function is (i? — T + ze)~^, and the scat- 
tering potential is Vop. As /c — > 0, the large-r behaviour 
of the first term on the RHS of Eq.([2|) gives the back- 
ground scattering length, and the second term yields the 



B. Energy spectrum in harmonic trapping 

Each particle, of mass M, moves in a harmonic po- 
tential (l/2)Mtj^r^. Making the usual transfomations 
to relative and CM co-ordinates, r = (ri ^ r2); and 
R = (ri -I- r2)/2, and their corresponding canonical mo- 
menta p, P, we obtain, for the noninteracting particles, 

^° = (2!^ + ^^^'"'^'^O + (C + ^-''") ' 

where Mem = 2Af, fi = M/2. Consider the rela- 
tive motion of these two trapped particles, interacting 
with a short range effective potential Vs{r). This is 
the single-channel equivalent of the potential W'^{r) of 
Eq.Q, with the energy dependence in the denominator 
being absorbed in the rapidly varying scattering length 
parametrized by the Zeeman splitting. The two-body 
Schrodinger equation in each partial wave / for the radial 
wave function ui(r) = rtpi^r) with energy Ei may be ex- 
pressed in dimensionless variables x — rj{\pl L), where 
L — ^h/{MLo), and r/; = 2Ei/{huj). It is given by 

d^ui m + i) 2Vs 2 

We may find the energy spectrum of the above equation 
without specifying the specific form of Vs by generalising 
a method first adopted by Jonsell @ for ^ = 0. Let the 
range of the short range potential Vs{r) be given by b. 
For r > b, taking Vg = 0. the solution of Eq.® is given 

by 



3 



ui = e 



, ,,1 f2l + 3-in 3 \ , , /-2/ + l-77( 1 

c\x'+'Mi — -J + ^;y] +c',x-'Mi -,^-i;y 



r 



(7) 



where y = x^, and M{a,j;z) is the confluent hyper- the wave function behaves as 
geometric function p^ . Since, for large z, M{a,j;z) 
behaves as 



rf7) 

M(a,7;z)~-P^z"-^e^ 
r(a) 



(8) 
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, r(3/2 + o 

-'ir((3 + 2^-770/4) 



r(i/2-0 



■r((l-2/-r;0/4) 

I 



^-(l+r„)/2ga;V2 
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for large x. In order to get a convergence solution, we 
have to have 



r(3/2 + or((i-2;-r?o/4) 
r(i/2-or((3 + 2;- 770/4) 



Since 



r(3/2 + 

r(i/2-0 



l\ [ (2^-1)!!]^ 
22^ 



(10) 



(11) 



we obtain 



^ = -(-1)' 



[(2^-l)!!]2r((l- 2^-7^0/4) 
22' r((3 + 2?- 770/4) 



(12) 

We now need to relate this ratio to the scattering length 
ai to determine the eigenvalue rji . Note that if the range b 
of Vs is very small, as we approach 7-^6+, the oscillator 
potential may be neglected. This assumption is crucial 
for our derivation, and is examined in some detail in 
Appendix A for I = 2. Thus ui{r) for positive energy 
may be regarded as the phase-shifted scattering solution 
due to Vs, given by 



ui{r) = Aikr[ji{kr) - ni{kr)t&n5i\ , 



(13) 



El = h^k'^ /M. For r > b, but stiU ~ 0, ui{r) behaves as 



ui(r) 
where 



B,r 



l + (2? + l)[(2/-l)!!]2 



tan 5i 



[kr) 



21+1 



Bi = 



Aik' 



(2;h-1)!! 



(14) 



(15) 



We match ui, given by Ea. (|T4)l with that obtained from 
Eq.([7]) for small x. The validity of this procedure for / = 2 
rests on the condition (|24p derived in Appendix A. For 



z Q, and 7 7^ 0, the Kummcr's functions M{a,^; z) is 
unity; hence the solution given by Eq.(I7]) behaves as 



ui{x) 



1 



'-2 -21-1 



When this is matched with Eq. p4|) . we obtain 
' tan 5i 



(2Z + 1)[(2/- 1)!!]2 



(V2fci) 



27+1 



(16) 



(17) 



Equating Eqsp^ [T7|) . and using the effective range ex- 
pansion ([1]). we obtain 



1 r((i-2/-7?o/4) 

V2r((3 + 2/-770/4) 



ai 



1 - ainrii/A 



(18) 



In the above, we have defined the dimensionless quanti- 
ties 

ai 



ai 



^27H 



n - L^'- V7 . 



(19) 



If we set I = 1 and L = Ir/V^ following Yip Q, we 
recover his result, that is, his Eq.(13); 

I'. 



-M2c)m = 8 



r((5-7n)/4) 

r((-l-7yi)/4) 



(20) 



where we have put ai = v and ri = 2c. Similarly, by 
setting Z = in Eg. p^ we recover the spectrum of the 
I = states as given by Busch Q . Note that Eq. has 
been obtained with no mention of any specific shape of 
the potential, and is valid for any short-range two-body 
potential. 



III. 



RESULTS AND DISCUSSION 



In Figs 1-3, we plot the energy spectrum given by 
Eq. as a function of the scaled scattering length 5/ for 
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1 = — 2, keeping the effective range f/ fixed. The same 
energy spectra look quite different when plotted against 
the inverse of the scattering length, l/o;. For 1 = 1, 
we reproduce Yip's result, and do not duplicate it here. 
Our Fig. 4 shows the energy levels for Z = 2, plotted as a 
function of 1/a;. Before discussing these spectra, we com- 
ment on the choice of the effective range parameter. For 
1 = 1, Yip had set the scaled effective range f i = -64/V2 
from experimental data . From a theoretical point of 
view, to see if this may be generated by a potential whose 
range is much smaller than the oscillator length L, we 
take a square-well potential. The shape of the potential 
should not matter unless the shape-dependent term in 
the effective range expansion (proportional to fc"*) affects 
the energy spectrum. We have verified that the spectra 
in Figs 1-4 remain virtually unchanged when this term 
is included. In the Appendix B, the analytical expres- 
sions for the scattering length and the effective range 
for any given partial wave are given. At a resonance, 
ai = ±oo, hence it follows from Eg. ([25)1 that j;_i(s) = 0. 
Therefore, from Eg. ([27]) we get vq = b, ri = —3/b, and 
r2 = -15/53. For Yip's choice of riL = n = -64/^(2), 
it follows that b/L = 6 ~ 1/15, fulfilling the condition 
that b << L. Unlike the case / = 1, we do not have 
guidance from experiment for the choice of r2, so we de- 
duce it from the square-well potential with 6 = 1/15, 
1/30, and 1/10. Furthermore, we find that for the 
square-well example, although the scattering length is 
highly sensitive to the choice of the strength parameter s 
defined by Eg. ([25)) . the effective range hardly changes 
as s is varied over a narrow range to accommodate 
the variation in the scattering length shown in Figs.(l- 



4). 
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FIG. 1: (Color online) Plot of the I = energy levels 
E in units of fkv versus the s-wave scattering length ao 
in units of the oscillator length L. The scaled effective 
range is fixed at 1/15. The plots change negligibly even 
for zero range. 
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FIG. 2: (Color online) The / = 1 energy levels vs the 
scaled p-wave scattering length a\. The dashed curves 
are for the effective range f\ = —30, while the continu- 
ous curves are for f\ = —45. 
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FIG. 3: (Color online) (a): The I = 2 energy levels vs 
the scaled d-wave scattering length. The dashed line is for 
f2 = —15 X 10'', while the continuous line is for r2 = —15*. 
For the excited states, the dashed and continuous lines can- 
not be distinguished on this scale. See text for the choice 
of the effective range parameters. The superimposed dots 
are the numerically calculated energies obtained by solving 
the eigenvalue equation (|6]) directly, (b): Amplified version 
for the excited states showing that the levels do not actually 
cross, but bend sharpely in Fig 3a. 



Figs. (1-3) show that away from the resonance (a; 
±(X)), the energy plots against di have similar shapes. 
But whereas the / = plots remain essentially unchanged 
when tq is varied over a wide range, the higher partial 



waves become more and more sensitive to the choice of 
the effective range. This is apparent from Figs. 2-3. In 
all three, the lowest energy state tends to —00 as di 0. 
When di ±00, the I ~ energy levels tend to the limit 
( 2?! 4-1/2), ?i = 0, 1,2.. in units of the oscillator spacing. 
By contrast, for Z 7^ 0, all but the lowest energy level go 
to the noninteracting values (2n -I- I + 3/2) in the zero 
range limit. All energies are in units of the oscillator 
spacing. At resonance, the lowest energy level for I ^ 
tends to zero as the range of the potential is decreased. 
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FIG. 4: (Color online) Plot of the excited state I = 2 energy 
levels vs the inverse of the scattering length in dimensionless 
units. This figure should be compared to Fig. 3 b. 



In Fig. 4, we see that the higher energy levels show 
strong bends away from the resonance, a feature even 
more pronounced than for the / = 1 case shown in Q. 

To confirm that the analytical result (|18p does produce 
the energy spectrum for / = 2 accurately for the range of 
the scattering parameters shown in Fig. 3(a), we solved 
Eq. (O numerically to determine 772 . For Vs , a square- well 
potential was taken with b = 1/15, and it was kept fixed. 
The scattering length 02 was varied as shown by chang- 
ing the depth of the potential. In Fig. 3(a), the dashed 
curves and lines are plotted using the analytical formula 
(fT5)) . The superposed dots show the results for the eigen- 
values obtained by the numerical integration of Eq.(l6l). 
The agreement is excellent, confirming the expectation of 
Appendix A. Table 1 shows, however, that the agreement 
with numerical integration becomes poor for the ground 
state as its energy approaches zero for large values of 0.2. 
This limitation is understandable from Eq. ((24|) derived 
in Appendix A and the discussion that follows it. 
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TABLE I: Comparison of energy 772 obtained numerically from Eq.© and from the analytical formula (|18|l . The range b) of 
the square-well potential is kept fixed at 1/15 



I = 2 


10'*a2 




-2 


-0.8 


-0.2 


0.2 


0.5 


1 


2 


100 


10^ 


10® 




E{hLo) 


Ground state 


N" 


0.1992 


0.4952 


1.9726 


-1.9772 


-0.7887 


-0.3945 


-0.1956 


-0.0022 


0.0017 


0.0017 




A' 


0.1975 


0.4939 


1.9762 


-1.9731 


-0.7894 


-0.3951 


-0.1976 


-0.0040 


-0.0001 


-0.0001 


First excited state 


N 


3.5002 


3.5010 


3.5012 


3.5001 


3.5001 


3.5001 


3.5002 


3.5002 


3.5002 


3.5002 


A 


3.5002 


3.5002 


3.5003 


3.5001 


3.5001 


3.5001 


3.5002 


3.5002 


3.5002 


3.5002 



"Numeric results 
''Analytical results 



111 summary, we have argued in this paper that it is 
legitimate to use the effective range formalism for higher 
partial waves in the presence of a coupling to a Fesh- 
bach resonance, even though these parameters may not 
exist for a long range 1/r^ interatomic potential. Next, 
we have derived an equation (|18p . valid under certain re- 
strictions for any partial wave, relating the eigenenergies 
to the effective range parameters. This is known to be 
be valid for / = and / = 1 as long as the range of the 
potential is much shorter than the oscillator length We 
have shown in this paper that it is also applicable for 
Z = 2, so long as the ground state energy is not too close 
to zero. Our Eq. fM]) . derived for the matching distance 
for I = 2 shows the limitation of the shape-independent 
parameters. At resonance, since both the terms on the 
RHS vanish, the analytical result becomes of limited 
validity. For I > 2, the restrictions are more severe. 

We are indebted to Dr.Takahiko Miyakawa for many 
useful discussions at Tokyo University of Science (TUS), 
where much of the work was done. We would also like to 
thank the referee for the incisive comments that resulted 
in our adding Appendix A, and a deeper understanding 



of the limitations of universality for I > 2. This work 
was supported by NSERC (Canada), and a grant from 
(TUS). 



IV. APPENDIX 

A. Appendix A 

The energy spectrum Eq.((T8]) resulted when the ratio 
Cj/c']^ from ([7]) with the oscillator potential present was 
equated to the corresponding ratio from the scattering 
solution without the confining potential. In doing 
so, we assumed that this matching was done at a small 
enough distance (still larger than the range b of Vs{r)) 
that the oscillator potential could be neglected. In this 
appendix, we justify this assumption by taking the spe- 
cific example oi I = 2. To do this, we need to write 
Eg. pop in more detail, and examine the neglected terms 
in the expansion of Eq.(I7]). For / = 2, for small x, this is 
given by 



J 



U2{x) 



m 5 



(21) 



It is in the coefficient 



of the x^ term that the To estimate the RHS, we use Eq.dT]), and put fc^ = 



oscillator presence is felt; without the oscillator, this co- 

2 

efficient is ^ . In Eq. (fTB|) , we neglect this term of order 
x^ , and yet retain the leading term x^ of the regular so- 
lution, even though x is small. This can only be justified 
if the matching distance x satisfies the inequality 



x^ > 



I- 

' 4 



(22) 



To check this, we substitute above the expression ([TT]) for 
c'2/c'i, which, for I = 2, gives the condition 



X > 



45 I tan (52 1 
T(7IfcI)5 



(23) 



We then get, for the condition 

1 1 

02 



'^E2 



(24) 



From Figs. 3 and 4. we see that i is of the order of 10'* 

(either sign), and f2 is —15 x 10'^ or ten times larger. The 
energy E2/hhj is of order unity. Unless there is some ac- 
cidental cancellation , the RHS of Eq. is very small, 
of the order of 10""*. On the LHS of Eq.dH]), x stands 
for the matching distance (in units of L, which is slightly 
larger than h = 1/15. Therefore the inequality condi- 
tion ([24|) is easily satisfied when the scattering length 
0,2 is small. This is not the case at resonance, however. 
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when 0,2 tends to infinity, and E2 for the ground state 
approaches zero. We then see that the RHS of Eg. ([24]) 
tends to infinity, and the inequahty condition cannot be 
satisfied. Our formula is no fonger accurate for the 
ground state in this situation. This is confirmed by our 
numerical calculation as shown in Table 1. 



In a given partial wave the expressions for the scatter- 
ing length ai and the effective range r; are given by 



B. Appendix B 



ai 



ji+i{s) 



(2;-l)!!(2/ + l)!!j,_i(s) 



(26) 



Consider an attractive square-well potential of depth 
Vb and range h. Define the strength parameter 



{^MV^/h? ) b 



(25) 



n 



(2/- 1)!!(2Z-H)! 

1,21-1 



21 - 1 



2; + lji-i(.s) 



ji+iis) 2l + 3\ji+i{s) 



3h (s) 



(27) 



[1] C. A. Regal, C. Ticknor, J. L. Bohn, and D. S. Jin, Phys. 

Rev. Lett. 90, 053201 (2003). 
[2] A. Marte et ai, Phys. Rev. Lett. 89, 283202 (2002). 
[3] S. Inouye et ai, Nature (London) 392, 151 (1998); Ph. 

Courteille et ai, Phys. Rev. Lett 81, 69 (1998). 
[4] T. Loftus, C. A. Regal, C. Ticknor, and D. S. Jin, Phys. 

Rev. Lett. 88, 173201 (2002). 
[5] T. Busch, B-G Englert, K. Rzazewski and M. Wilkens, 

Foundations of Physics, 28, 549-559 (1998). 
[6] S. Jonsell, Few-body Systems 31, 255-260 (2002). 
[7] P. Shea, B. P. van Zyl, and R. K. Bhaduri, Am.J.Phys. 

(in press) 2009. 
[8] T. Stoferle et al, Phys. Rev. Lett. 96, 030401 (2006). 
[9] S.-K. Yip, Phys. Rev. A78, 013612 (2008). 
[10] C. Ticknor, C. A. Regal, D. S. Jin, and J. L. Bohn, Phys. 

Rev. A69, 042712 (2004). 
[11] S.-K. Yip, arXiv:0801.0636y l (unpublished). 



[12] N. F. Mott and H. S. W. Massey, The Theory of Atomic 
Collisions (Oxford University Press, Third Edition, 1965) 
p.50. 

[13] L. D. Landau and E. M. Lifshitz, Quantum Mechan- 
ics, (Addison- Wesley, Reading, Massachusetts) (1958), 
p.405. 

[14] Bo Gao, Phys. Rev. A58, 4222 (1998). 

[15] T. Kohler, K. Goral and P. S. Julienne, Rev. Mod. Phys. 
78 1311 (2008) . 

[16] M. Block and M. Hohhaus, Phys. Rev. A65 052102-1 
052102-4 (2002). 

[17] C. Cohen-Tannoudji, in Lectures on Quan- 
tum gases, Institut Henri Poincare, Paris, 
April 2007 (unpublished). PDF available at 
|http: //www. phys. ens.fr/~castin/progtot. htm l 

[18] B. R. Levy and J. B. Keller, J. Math. Phys. 4, 54 (1963). 



